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The collective strength of a system of fibers, each having a failure threshold drawn randomly from a 
distribution, indicates the maximum load carrying capacity of different disordered systems ranging 
from disordered solids, power-grid networks, to traffic in a parallel system of roads. In many of 
the cases where the redistribution of load following a local failure can be controlled, it is a natural 
requirement to find the most efficient redistribution scheme, i.e., following which system can carry 
the maximum load. We address the question here and find that the answer depends on the mode 
of loading. We analytically find the maximum strength and corresponding redistribution schemes 
for sudden and quasi static loading. The associated phase transition from partial to total failure 
by increasing the load has been studied. The universality class is found to be dependent on the 
redistribution mechanism. 


The response of an ensemble of elements having ran¬ 
dom failure thresholds plays a crucial role in the break¬ 
down properties of heterogeneous systems [l[ . Such sys¬ 
tems can be, for example, disordered solids under stress, 
power-grid networks carrying currents, roads carrying car 
traffic, or redundant computer circuitry etc. The fail¬ 
ure properties of these systems are strongly dependent 
upon the load redistribution mechanism following a local 
breakdown. While in some cases (e.g. solids under stress) 
such mechanisms are properties inherent to the system, 
in many other cases (e.g. power-grids Q, traffic controls 
i) it is a matter of design that can be optimized so as 
to achieve the most robust configuration (see e.g. i). 
Such robustness properties of networks connecting ele¬ 
ments with varying failure thresholds under targeted or 
random attacks have received substantial attention 

In this Letter, we address the question of maximiz¬ 
ing the strength of a disordered system by finding the 
most effective redistribution scheme analytically using 
the fiber bundle model as a generic example. This model 
was introduced Q, indeed, to estimate the strength of 
cotton in textile engineering. Since then, it has been ex¬ 
tensively studied in the context of distribution of failure 
strength and failure time of disordered materials under 
tensile loading or twist by viewing fibers as elements of 
the disordered solids having a finite failure threshold (or 
even a finite lifetime dependent of loading) [§-[^ (see 
[ 2 ^ for a review). 


Conventionally, the fiber bundle model is viewed as 
a set of parallel fibers, having failure thresholds ran¬ 
domly drawn from a distribution (say, uniform in [0 : 1]), 
clamped between two horizontal plates. When the bot¬ 
tom plate is loaded, some of the weak fibers break, and 
their load is redistributed among the surviving fibers, 
which may in turn break or survive depending on their 
failure thresholds. The load transfer rule, generally a 
function of distance from the broken fiber, plays a vital 
role (see e.g., [HIH). While substantial attention has 


been concentrated in revealing the transition properties 
and related dynamics at and near the failure point, much 
less is known about tuning the redistribution scheme to 
maximize the strength of the system. In many applica¬ 
tions, as mentioned above, this is a matter of design and 
hence remains an important question. 


We consider the system as simply a set of intercon¬ 
nected elements and the load transfer is heterogeneous 
and dynamic in the sense that it may depend on the 
threshold and also on the instantaneous values of the 
loads on different elements. Over the years, examples of 
time dependent failure, which could sometimes be trans¬ 
lated to some sort of threshold dependence, were studied 
(see for example 1^, 18, 2^). However, it is hard to make 
explicit comparisons here. 


Before entering into the details of the heterogeneous 
scheme, we note that the loading process can be broadly 
divided in two ways: sudden loading and quasi-static or 
gradual loading. We also note that the external loading 
is always uniform although redistributions can lead to 
inhomogeneities in load carried by individual fibers. 

In the case of sudden loading, a load W is applied uni¬ 
formly on the initially intact and unloaded system in a 
step function, consequently the load per fiber initially is 
a — W/N {N being system size, i.e. number of fibers). 
The redistribution schemes then follow the initial break¬ 
ing, leading to a stable state which is either partially 
damaged or has completely failed. Irrespective of the 
redistribution scheme followed, an upper bound of fail¬ 
ure threshold (say, dm) can readily be estimated for this 
case. On application of stress dm to each fiber, those hav¬ 
ing thresholds < dm will immediately break. The best 
possible scenario is that no more fibers break and they 
share the total load and carry load to their highest ca¬ 
pacity. The average threshold of the remaining fibers is 
(1 + cr„i)/2 and the fraction of surviving fibers is (1 — dm) 
(threshold distribution uniform in [0 : 1]). The condition 
of maximum strength becomes (1 — dm ) (1 + o'm ) / 2 = dm 
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i.e., am, = \/2 ~ 1- For arbitrary threshold distribution, 

OO 

the condition is / fp{f)df = am, where p{f) is a thresh- 

<7m 

old distribution function, i.e. total remaining threshold 
should be equal to total load (see Supplemental Mate¬ 
rial). This limiting value is obtained without assuming 
any particular load redistribution scheme and it is worth¬ 
while to find out at least one scheme which can lead to 
this limiting behavior. 

The uniform load redistribution is a homogeneous 
scheme which can give a loading capacity equal to 1/4 
[ 2 ^, much less than the maximum value derived above. 
Hence we use a heterogeneous scheme to explore the pos¬ 
sibility of obtaining a failure threshold larger than 1/4. 

Particularly, we have considered the case where the 
load received by the i-th fiber in a given redistribution 
process is proportional to {fi — aif’, where ai and fi 
are respectively the load and failure threshold of the f-th 
fiber. The limit 5 = 0 is the well studied 23| uniform load 
redistribution. We have numerically obtained the critical 
loads for different b values, which also depend on whether 
the load is applied suddenly or quasi-statically. We could 
analytically calculate the limiting cases of critical load 
and have found that the most robust condition (highest 
value of critical load) arises for 6 = 1 in sudden loading 
and 5 —7> OO for quasi-static loading (effectively for 5 ~ 10 
or above). 

Furthermore, from the point of view of phase transi¬ 
tion between partial and complete failures, we have found 
that for 0 < 6 < 1, the order parameter exponent value 
remains the same. But for 5 > 1 it is different, showing a 
different universality class than that of the uniform load 
sharing one. The dynamics, and hence the system size 
dependence of the relaxation time near the critical point 
depends strongly on the loading mechanism. 

A natural expectation for a better scheme is to assume 
that the stronger fibers should get a higher share. How¬ 
ever, the effectively stronger fibers are those which have 
higher difference between the failure threshold and the 
load it presently carries. Hence, we consider that the 
load received by the i-th fiber having stress ai and fail¬ 
ure threshold fi is Xi = A{fi — aifiUe, where We is the 
excess load to be redistributed in a given step and 6 is a 
parameter. In general, the redistribution can be a many 
step process, but let us consider the first step only, such 
that ai = a and We = a'^N (since in the first step aN 
number of fibers will break each having load a). The 
normalization (or, load conservation) condition reads 


^ J [A{ff - afa^N + u] P{ff)dff = aN, 

i^survivingfibers ^ 

( 1 ) 

with aN being the applied load and P{f) oc 1/(1 — a) 
(this is essentially the normalized distribution of thresh¬ 
olds of remaining fibers), giving A = (6-f l)/(Af(l—cr)^+^) 



FIG. 1: (Color online) The phase diagram in the b — ac 
plane is shown for (a) sudden and (b) quasistatic 
loading for various fractional errors in the knowledge of 
the threshold values of the individual fibers (curves are 
from top to below for e = 0.0,0.1, 0.2, 0.3, 0.4, 0.5, 0.75). 
The upper bounds for both cases are shown which are 
reached for 6 = 1 (a) and 6 —^ 00 (b). 


(see Supplemental Material for calculation). Hence the 
load received by the f-th fiber in the first redistribution 
step, having failure threshold fi following the sudden ap¬ 
plication of load per fiber value a is 

( 2 ) 


Now, the redistribution will occur only once (no recursive 
dynamics) if after the first redistribution all fibers have 
force lower than their respective thresholds i.e., x{fi) + 
a < fi for all i £ surviving fibers, (i) For 6=1, this 
inequality leads to 




■ 2^2 



> a 


■ 2^2 



(3) 


which is generally satisfied (since for surviving fibers 
fi > a) upto a critical CTc where 2a\/(\ — affff = 1 or, 
ac = — 1- But this is the maximum possible sus¬ 

tainable load per fiber value, implying for 6=1, the 
maximum strength is achieved and that is a single step 
redistribution (see Fig. [!}; (ii) For 6 > 1, the condition 
for recursive dynamics leads to a lower bound for fail¬ 
ure strength, given by ai = (see Supplemental 

Material for derivation). 

The exponent values of the failure transition can be 
estimated for different 6 values. Particularly, the order 
parameter exponent value can be analytically found for 
6 = 1 and is known for 6 = 0. The well studied case of 
6 = 0 gives the order parameter variation as O = U{a) — 
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FIG. 2: (Color online) The variations of the derivative 
of order parameter are shown for different b values for 
(a) sudden and (b) quasi-static loading. The quantity is 
inverse of 1 — C (see Eq. (|H)), which is also used as an 
alternative order parameter. These indicate that the 
order parameter exponent value /3 = 1/2, for & < 1. For 
b> 1, /3' = 0, i.e. the curves gradually become 
horizontal. 



b 


FIG. 3: (Color online) The main panel shows the 
variations of maximum relaxation times (pre-critical 
and post-critical) with b for different system sizes. As 
6 —)► 0, the behavior is maintained (see inset). But 

as b is increased, the power-law behavior is lost. For 
b > 1, system size dependence is almost non-existent. 
These data are for sudden loading. For gradual loading 
similar behavior is seen, as the two versions are 
equivalent as far as the dynamics on or around critical 
point is concerned. 


U{ac) = C{ac — cr)^, with (3 = 1/2 [^, where U{a) 
is the fraction of surviving fibers under a load per fiber 
value cr. This definition is, however, difficult to follow 
in the numerical measures where the quantity U{(Jc) is 
not known exactly. A different order parameter measure, 
called the branching ratio C, was proposed in (see also 


( 2 ^ 1 . It was defined as C = 1 — n^/ ^ n*, where ng is 

t=o 

the number of failed fibers immediately following a load 
increase (necessarily small amount) and the sum is the 
total number of fibers failed due to these initial failures. 
Then the quantity 1 — C near the critical point follows 
1 — ^ ^ (ctc — <jY ■ The exponents [3 and (3' need not be 

tmax 

equal in general. The quantity ^ rit is essentially the 

t=o 

change in the number of fibers following a small increase 
in the load, or in other words is proportional to \dU/da\. 
But near the critical point 

- 1 - c-- (CTc - cr)^"\ (4) 


dU 

da 


This leads to the relation (3 = 1 — (3'. Incidentally, for 
b = 0, (3 = 1/2, giving (3 = (3' in that case. However, 
for 6 = 1 the redistribution process has been shown to 
be single step and U{a) = 1 — a. Therefore, O = U{a) — 
U(ac) = {oc — a), giving (3 = 1. On the other hand, since 
there is no avalanche in this case until the failure point, 
1 — C = 1 for a < Uc and 0 at cr = CTc (one may think this 
as a zero value of the exponent (3'). 

Interestingly, the exponent (3 seems to have a value 
close to 0.5 for 6 < 1 (see Fig. [5] (a); (b) shows the same 
for gradual loading). For 6 > 1, /3 remains equal to 1, 
however, the value of \dU/da\ shows a sudden change at 
CT;, where the redistribution process starts extending to 
more than one step. 

Another way to probe the critical behavior is to look at 
the relaxation time at (or very near to) the failure point. 
It is known that at the critical point, the relaxation time 
scales as T ^ A^“, with a = l/3for6 = 0 25, 2^. But 
for b > 0 the power law scaling seems to be lost, with 
r becoming almost independent of for 6 > 1 (see Fig. 
IS. Here we have shown the behavior of both the pre- 
critical relaxation time (the relaxation time in the last 
stable configuration of the model) and the post-critical 
relaxation time (the relaxation time in the first condition 
of complete failure). 

Apart from the sudden loading limit discussed so far, 
the other extreme is the quasi static loading. Here at 
every step the load is uniformly increased on all the sur¬ 
viving fibers until the weakest one breaks. The loading 
is then stopped until the system goes to the next stable 
state by local failures and subsequent stress redistribu¬ 
tion. After the relaxation to a stable state the system is 
further loaded to continue the dynamics (system relax¬ 
ation assumed to be much faster than external loading 
rate). By observation it is seen that the failure threshold 
increases with b. For a very large b, any redistribution 
of the load will essentially be eaten-up by the fiber(s) 
with largest {fi — at) value(s). Therefore, in the situa¬ 
tion where the load is gradually increased upto the point 
where a fiber with strength / has broken, the load pro¬ 
file will be like the one schematically depicted in Fig. 0] 
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FIG. 4: (Color online) (a) shows the schematic load 
profile in 6 —>■ 00 limit when the gradual loading is done 
upto some threshold value. The fibers are arranged in 
increasing order of their threshold value. The black 
triangle on the left are the broken fibers. The grey area 
denotes the load. There is a uniform external loading. 
The redistributed load of the broken fiber is dumped on 
the right hand corner triangle, because for those fibers 
fi — ai were highest, (b) shows the extreme limit, 
beyond which all fibers will break together. Simulation 
results for 6 = 10 supporting the above picture is shown 

at (c). 


(a). The redistributed load will be piled up in the right 
hand corner in the form of a triangle and there will be a 
uniform loading for the other fibers, which is due to the 
uniform load increase mechanism. Now, the extreme end 
of this process is the one depicted in Fig. 0] (b). Here 
/m = 1/2 and beyond this point all remaining fiber has 
(fi — ai) values equal to each other. Hence, a further 
increase in the load will result in the complete collapse 
of the system. This is, therefore, is the point where the 
system carries maximum load using gradual loading. The 

total load (per fiber) carried by the system at this point 
1 

is f a da = 3/8. Therefore, maximum failure threshold 

1/2 

for gradual loading is 0 '/’’“'^ = 3/8, for uniform thresh¬ 
old distribution in [0 : 1], and is obtained in the limit 
& —^ 00 . In practice, for 6 ~ 10 the maximum value is 
almost achieved. Fig. 0] (c) depicts the situation when 
the simulation is done for 6 = 10. 

The universality class in this case is determined as be¬ 
fore. The derivative of the order parameter diverges with 
|(Jc — cr| with an exponent value 1/2 as long as & < 1 
(see Fig. 0]upper panel). In this case, unlike the sudden¬ 
loading one, 6=1 too seems to belong to this universality 
class, at least within the numerical accuracies. For 6 > 1, 
the divergence gradually flattens as the critical point is 
approached, indicating a change in the universality class. 


as before. The behavior of the relaxation time remains 
the same (not shown) as in the case of sudden loading. 

The maximum achievable strength in sudden loading 
is obtained for 6 = 1 following the load redistribution 
method proposed in this Letter; however this scheme re¬ 
quires the detailed knowledge of the failure thresholds of 
each element. In real situations for controlled processes 
this knowledge can be rather imprecise. The question 
is how effective these load redistribution strategies are 
in those cases? Let us say the imprecisely known value 
of the threshold of a fiber /“ can be anywhere between 
fi{l — e) < /“ < fi{l + e) with uniform probability, where 
e is a constant (< 1) and fi is the actual threshold (our 
results so far are for e = 0). We set the additional rule 
that the fibers already carrying loads equal to or higher 
than their assumed thresholds will not get any further 
share of load from redistribution, which can happen when 
the assumed value is less than the actual value (due to 
this rule, even for 6 = 0 the critical load is higher than 
1/4, since, even for an imprecise knowledge, generally the 
weaker fibers do not get load in redistribution). On the 
other hand, if the assumed value is larger, it will break 
when the actual value is exceeded. 

It turns out, even for e > 0, this strategy can give 
better result than uniform load redistribution for both 
sudden and quasistatic loading (see Fig. [T]) with CTc signif¬ 
icantly higher for a considerable range of value of 6 > 1; 
the range decreasing with higher values of e. For the 
quasi static loading case, the emergence of a shallow peak 
in CTc occurring at a finite value of b is to be noted. These 
observations indicate the effectiveness of the present re¬ 
distribution scheme in enhancing the upper critical value 
of the load even with partial knowledge. 

While we have noted marked difference between the 
two extreme loading mechanism, there are instances for 
time dependent fiber bundles [l9|, [22 1, where these two 
limits can tend to converge. Detailed discussions on these 
are beyond the scope of the present manuscript and will 
be taken up elsewhere [s^ . 

In conclusion, the threshold dependent load redistribu¬ 
tion scheme proposed here can lead to significantly higher 
strength in the fiber bundle model, compared to the uni¬ 
form load redistribution mechanism usually followed. For 
threshold distribution uniform in [0 : 1], the maximum 
strength (load per fiber value prior to catastrophic fail¬ 
ure) for sudden loading is found to be \/2 — 1, which 
can be achieved for load redistribution proportional to 
{fi — ai)^, with 6=1. For quasi static loading, the 
maximum strength is 3/8 and is achieved in the limit 
6 —>■ c». While for sudden loading the maximum achiev¬ 
able strength is obtained using this scheme, for qua¬ 
sistatic loading whether there exists a more efficient redis¬ 
tribution scheme that can lead to a higher value of max¬ 
imum strength, remains an open question. We also find 
that the universality class is different from the homoge¬ 
neous loading case for 6 > 1 for both the loading schemes. 
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The results remain qualitatively valid for other threshold 
distributions (Weibull, Gumbel, Gaussian) and are also 
useful when the individual thresholds are known only ap¬ 
proximately. The precise bounds obtained here and the 
corresponding load redistribution mechanisms can help 
in maximizing the strengths of complex interconnected 
systems with an ensemble of elements each having finite 
failure threshold. Of course, this opens up the challeng¬ 
ing question of practical implementation, which will de¬ 
pend on the specific details of the system considered. 
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Maximizing the strength of fiber bundles 
under uniform loading: supplemental 
material 


DERIVATION OF LOWER BOUND OF 
CRITICAL LOAD FOR SUDDEN LOADING 


One can show that for & > 1, if there is a recur¬ 
sive dynamics that continues beyond the first redistribu¬ 
tion step, then it initiates from the fibers having highest 
thresholds. As mentioned in the main text, the incre¬ 
ment of the load on a fiber with failure threshold fi in 
the first redistribution step, following a sudden applica¬ 
tion of load tr is 

*(/<) = ® 

Hence the total load on a fiber with threshold fi becomes 
CT + x{fi). Since all fibers having threshold on or below 
a have already broken immediately after the application 
of load, all remaining fibers are of threshold higher than 
a. Hence we can write 


fi = a Ci. 


( 6 ) 


First we show that for very small e^, the fiber survives. 
The condition for survival for a fiber is (dropping the 
index i) 


x{f)+a < f 


(7^(6 -bl) f, 

^ -e'’-b (7 < cr-b e 


(I - a)&+i 

cr^(&-b 1) j,_i 

(I-a)'>+i 


< I. 


(7) 
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FIG. 5: (a) The critical stress for sudden loading is 
plotted (this corresponds to the e = 0 line in Fig. 1 (a)) 
along with the lower-bound (Eq. |9]) for 6 > 1. The lower 
bound essentially marks the onset of recursive dynamics 
and hence matches with the critical load for 6 = 1 
where the process is single step, (b) The stress profiles 
on the fibers (arranged in ascending order of their 
failure thresholds) after the first redistribution step for 
different a values and for 6 = 2. For a > ai (~ 0.366 in 
this case), the fibers with high threshold values will 
break, indication the onset of recursive dynamics. 


This gives, 




y/TTb-1 

b 


(9) 


which is essentially a lower bound of CTc for 6 > 1. These 
are also supported by numerical simulations (see Fig. [S] 
(a)). An estimate of this lower bound may be useful for 
safe designs. 


CALCULATION OF THE NORMALIZATION 
CONSTANT A 

The normalization condition (Eq. (1) in the main text) 
reads 

1 

^ f [A{f, - afa^N + a] P{f,)df, = <jN, (10) 

while the sum is performed on the surviving fibers (i.e. 
i G S) following an application of load a on each fiber 
(hence the limit in the integral), the right hand side is the 
total applied load, which is to be conserved. Performing 
the integral 

E 

ies 

( 11 ) 

where we have used P{f) oc 1/(1 — cr) for the remaining 
fibers. Now the sum gives a factor (1 — a)N, which is the 
number of remaining fibers. Putting that in the above 
equation we have 


A 




6 + 1 


a'^N + (j{\ — a) 


1 — a 


= (tN, 


^ (b+1) 

fV(l - a)'’+i ’ 

which is used in the main text. 


( 12 ) 


This is satisfied for small e and when 6 > 1 (and a not 
close to unity, which would in any case break the entire 
system). Hence the breaking of the fibers will begin from 
the highest threshold value (unity in this case). This 
picture is also supported in Fig. [5] (b), where we have 
plotted the load profiles on the fibers (arranged in as¬ 
cending order of their threshold values) after the first 
redistribution step for different tr values. As can be seen 
from the figure, there is a value ct;, on or above which 
more fibers will break in this step, marking the onset of 
recursive dynamics. The value of ct; can be estimated 
from the condition that it is the load for which the fiber 
with highest threshold value (unity) will break after first 
redistribution. Hence 

( 8 ) 


CALCULATION OF UPPER BOUND OF 
FAILURE THRESHOLD FOR SUDDEN 
LOADING 

In the main text, the upper bound of failure threshold 
(am) was derived for the threshold distribution uniform 
in [0 : 1]. Here we provide the general expression for 
arbitrary threshold distribution and then focus on the 
specific type of uniform distribution (between [a : 1 — a]) 
for which the mean is fixed at 1/2. The a —>■ 0 limit is the 
one we studied in the main text, am has non-monotonic 
variation with a. 

First consider an arbitrary threshold distribution func¬ 
tion p(f) defined in [0 : oo]. If a load am is applied uni¬ 
formly on all fibers, then all fibers having threshold val¬ 
ues less than am will break immediately. Now, the best 
possible scenario (with regard to higher failure threshold) 


x(fi = 1) + ai = 1. 
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FIG. 6: (a) The uniform threshold distribution in the 
range [a : 1 — a] is shown with a schematic diagram. 
The maximum critical load am is also indicated, (b) 
The variation of am with the parameter a is shown. 
The a —>■ 0 limit is discussed in the main text. The 
maximum value occurs at a —>■ 1/2 limit, which is the 
limit of vanishing disorder (delta function limit). The 
non-monotonicity arises out of the competition of less 
number of broken fiber on application of load am and 
the presence of smaller number of stronger fiber as a is 
increased. 


is the one where no further fiber breaks and all of them 
carry load upto their maximum capacity. This means 


J fpif)df = am- (13) 


Solution of this equation will yield the maximum failure 
threshold for sudden loading am for an arbitrary thresh¬ 
old distribution. 

Now let us focus on the uniform threshold distribution 
with its mean fixed at 1/2 (see Fig. [n](a)). This is a gen¬ 
eralization of the particular distribution we considered in 
the text. The distribution has the limits [a : 1 — a]. So, 
the normalized distribution is 


Pif) = 


Applying Eq. (ITTTl) we have 


1 

1 - 2a' 


(14) 


l-g-crml-g + o'm 
1 - 2a 2 


— 1 


(15) 


of course the upper limit of the integral here is 1 — a. 
Solving it one has 

am = 'J 5a^ — 6a -I- 2 -|- 2a — 1. (16) 

The a —^ 0 limit gives fTm(O) = \/2 — 1, as discussed in 
the main text. The plot of am versus a is shown in Fig. 
[S] (b). In the a —>■ 1/2 limit, the threshold distribution 
becomes a delta function centred at 1/2. That is of course 
the highest value of am ■ But before that it shows a non¬ 
monotonicity with a. This is because, increasing a brings 
a competition between less number of broken fibers with 
less number of strong fibers present. 














